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We establish an explicit correspondence between ergodicity breaking in a system described by 
power-law tail distributions and the divergence of the moments of these distributions. 
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Statistical mechanics is a combination of the law of 
large numbers and the laws of mechanics. Since its foun- 
dation in the late 19th century, this theory has been ex- 
tremely successful in describing equilibrium and nonequi- 
librium properties of a very large number of macroscopic 
systems p|. In the last decade, however, a new class 
of systems that do not obey the law of large numbers 
has emerged 0, El 01 ■ The behavior of these systems is 
dominated by large and rare fluctuations that are charac- 
terized by broad distributions with power-law tails. The 
hallmark of these statistical distributions, commonly re- 
ferred to as Levy statistics p], is the divergence of their 
first and/or second moment. 

The question we address in this paper is how ergodicity 
is affected in systems described by power-law tail distri- 
butions with diverging moments. Ergodicity is a cen- 
tral concept in statistical physics and is usually stated 
by saying that ensemble average and time average of 
observables are equal in the infinite-time limit [g, Q- 
The ergodic hypothesis has recently been investigated ex- 
perimentally in two different systems governed by Levy 
statistics in time: fluorescence intermittency of nanocrys- 
tal quantum dots Q and subrecoil laser cooling of atoms 
[3- Both experiments have found that Levy statistics 
induces ergodicity breaking. A precise understanding 
of the connection between divergent moments and non- 
ergodicity is thus of high interest. A common feature 
of the systems in the above experiments is their non- 
stationarity. This in in contrast to the system we propose 
to study, namely atomic transport in an optical lattice, 
where a steady state does exist. Nonetheless, we will 
show that this system can exhibit nonergodic behavior. 

An optical lattice is a standing wave light field formed 
at the intersection of two or more laser beams [l0[ . Due 
to the spatial periodicity of the potential, an optical lat- 
tice is similar in many respects to a solid state crystal. 
The main advantage of an optical lattice, however, is its 
high tunability: both the period and the amplitude of 
the optical potential can be modified in a controlled way. 
This has the interesting consequence that the exponents 
of the power-law distributions appearing in this system 
are not fixed, as in most systems, but can be varied con- 
tinuously, allowing the exploration of different regimes. 

The motion of atoms in a one-dimensional optical lat- 
tice formed by two counterpropagating laser beams with 
linear perpendicular polarization can be described, after 
spatial averaging, by a Fokker-Planck equation for the 



semiclassical Wiener function W(p,t) 013113, 
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The momentum-dependent drift and diffusion coeffi- 
cients are respectively given by, 



K(p) = - 



ap 



1 + ip/Pc 



D(p) = D a 



1 + ip/Pcf 
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The drift K(p) corresponds to a cooling force with fric- 
tion coefficient a, while the diffusion coefficient D(p) rep- 
resents stochastic momentum fluctuations and describes 
heating processes. It is worth noticing that for large mo- 
mentum, the range of the drift is limited by the capture 
momentum p Cl while the range of the fluctuations is not. 
The kinetic equation Q is valid in a regime where (i) 
the atomic momentum is large, p 3> Tik, where k is the 
wave number of the laser field (this defines the semiclas- 
sical limit), (ii) the saturation parameter is small, s <C 1 
(this corresponds to low laser intensity) and (iii) the ki- 
netic energy of the atoms is large, p 2 /2m 3> Uq, where 
Uq is the depth of the optical potential (this last condi- 
tion allows spatial averaging) . The stationary solution of 
the atomic transport equation Q which satisfies natural 
boundary conditions, W s (p — > ±oo) — + 0, is of the form, 
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where Z is a normalization constant. The momen- 
tum distribution J3J has an asymptotic power-law tail, 
W s (p) ~ \p\~( a Pc)/ D o ( w ith an exponent that can be ex- 
pressed in terms of the potential depth Uq and the recoil 
energy E R as (ap 2 c )/D = U /(22E R ) [ll|. The statis- 
tical properties of the distribution J3J can therefore be 
easily changed from normal statistics for Uo > 66-Er, to 
Levy statistics for Uo < 66-Er, by simply modifying the 
depth of the optical potential. In particular, one should 
note that the second moment, (p 2 ) = J dpp 2 W s (p), be- 
comes infinite when Uq <66Er. In this regime the mean 
kinetic energy of the system, Ek — (p 2 )/2m, diverges, 
clearly signaling unusual thermodynamic behavior. 

Transport in an optical lattice has a number of attrac- 
tive features that make it an ideal case study of the ther- 
modynamical properties of systems described by power- 
law distributions. On the one hand, the atomic trans- 
port equation Q has been derived from the microscopic 
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Hamiltonian that describes the interaction with the laser 
fields and the quantities that appear in this equation can 
be expressed in terms of the microscopic parameters of 
the quantum-optical problem [TT| . Moreover, the regime 
defined by conditions (i) to (iii) has been implemented 
experimentally and the divergence of the kinetic energy 
below a given potential threshold has been observed [l4| . 
Finally, Eq. is an ordinary linear Fokker-Planck equa- 
tion, meaning that standard methods of stochastic cal- 
culus can be used to analyse the problem. This is in 
contrast to most systems with power-law distributions 
that are often described by nonlinear or fractional ki- 
netic equations |l5j |. In particular, for the case of van- 
ishing Di (which we shall consider in the sequel 
Eq. JQ) corresponds to a random process driven by ad- 
ditive Gaussian white noise. The fact that power-law 
fluctuations with infinite variance occur here in a sys- 
tem subjected solely to additive noise is worth empha- 
sizing. The physical mechanism which gives rise to di- 
vergent fluctuations in systems with multiplicative noise, 
where the noise strength is proportional to the stochas- 
tic variable, is well-known [T3: it is based on a positive 
feedback loop that leads to the amplification of the noise 
as the value of the random variable increases. On the 
other hand, in the present situation, where the noise in- 
tensity is independent of the random variable, a different 
positive feedback mechanism is at work, based on the 
steady decrease to zero of the friction force as the value 
of the momentum increases. The appearance of infinite 
momentum fluctuations in this system is thus a striking 
illustration of the complex behavior that can result from 
the subtle interplay of the noise and the nonlinearity of 
the drift. In the following, we establish a correspondence 
between the divergent moments of the power-law distri- 
butions of the system and nonergodic behavior. 

We begin by transforming the Fokker-Planck equation 
(|T|) into a Schrodinger-like equation by writing W(p, t) — 
Wsijp) 1 /" 2 x ip(Pyt) [iSl - The function ?p(p,t) satisfies the 
imaginary-time Schrodinger equation, 



dtfj 
~dt 



-D 



dp 2 



V(p)ip = Hip 



(4) 



with the potential V(p) = K'(p)/2 + K(p) 2 /(4D ). For 
the drift coefficient J2J), this potential reads 



V(p) = 



ap 2 c p 2 (ap 2 c + 2D )-2D Q pl 
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This transformation reveals the fundamental difference 
between the case of finite p c and the case of infinite p c , 
where Eq. JQ) reduces to the familiar Ornstein-Uhlenbeck 
process with linear drift. For infinite p c , the Schrodinger 
potential V(p) asymptotically increases asp 2 , whereas for 
finite Pc, it asymptotically decreases as 1/p 2 . By contrast, 
the Fokker-Planck potential, $(p) = — J Q P dp' K(p') 7 
is confining for any value of the capture momentum. 
As a consequence, the spectrum of the Hamiltonian, 
Hipk(p) = Ekipk(p), is discrete in the former case, while it 



is continuous, except for the discrete ground state ipo(p), 
in the latter. In both cases, the stationary momentum 
distribution is given by the square of the ground state 
eigenf unction, W s {p) — ipoip) 2 - Interestingly, we note 
that a real-time Schrodinger equation with a potential of 
the form © has recently been considered in Ref. . To 
simplify the discussion, we now adopt rescaled variables 
for which a = p c — 1 and -Do = D; the noise intensity D 
being then the only remaining parameter in the problem. 
It will also be convenient to divide momentum space into 
a low-momentum region \p\ < 1, where the drift is ap- 
proximately linear, K\(p) ~ —p, and a high-momentum 
region \p\ > 1, with a drift K^ip) — —1/p- The anoma- 
lous dynamics of the system is completely determined by 
the high-momentum region. The eigenvalues and eigen- 
functions of the Hamiltonian H in this region are given by 
E k = Dk 2 and ipk(p) = y/p[ciJv(kp) + c 2 Y v (kp)], where 
J u (p) and Y v (p) are the Bessel functions of the first and 
second kind of order v = (D + 1)/2D. The constants C\ t2 
are fixed by the boundary conditions. 

When discussing the ergodicity of a system, one is typ- 
ically not interested in the trajectory in the full space- 
space, but often in the projection of the trajectory onto 
some subspace of relevant variables Q], in the present 
case momentum. A criterion for the equality of ensemble 
average and time average of a dynamical quantity A is 
then provided by the condition 



Am 



((a -(A)) 2 ) 



when t 



00 



(G) 



Here A = t 1 f Q dr A(p(r)) is the time average of the 
observable A and (A) = J dp A(p)W(p,t) denotes its 
ensemble average. In the infinite-time limit, the lat- 
ter tends to the stationary ensemble average (A) s — 
J dp A(p)W s (p). A system that obeys © is said to be 
ergodic in the mean square sense. In order to determine 
the ergodicity of our system, we thus need to compute 
the long-time behavior of the covariance, 



t pt 
dti dt 2 



(A(p( tl ))A(p(t 2 ))) 
(AWh))) (A(p(t 2 ))) . (7) 



This can be done by applying the usual theory of stochas- 
tic processes ^^|. The two-time correlation function 
(A(p(ti))A(p(t2))) is defined by the integral 

(A(Kii))A(p(i 2 ))) = 

dp 1 dp 2 A(pl)A(p 2 )W 2 (pi,t 1 -p 2 ,t 2 ) ,(8) 



where W 2 (pi 1 ti;p 2 ,t 2 ) is the two-point joint probabil- 
ity density function. Since the Fokker-Planck equation 
describes a stationary Markov process (for any value 
of p c ), the probability distribution W 2 {pi,ti,p 2 ,t 2 ) de- 
pends only on the time difference \t\ — t 2 \ and can be 
expressed in terms of the eigenvalues and eigenfunctions 
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of the Schrodinger equation in the form, 
W 2 (p 1 ,t 1 ;p 2 ,t 2 ) = ip (pi)ip (p 2 ) x 
(MPi)MP2) + J™ dk MPi)MP2) e- E ^-^\) .(9) 

Combining Eqs. @, J8J and J^J and introducing the vari- 
able t — t 2 — t\, we arrive at 



2 f* 



(t-r)C A (r) 



where the function CU(r) is given by 



Ca(t) 



dk 



dpA(p) ipo(p)ipk(p) 



-E k r 



(10) 



(11) 



The infinite-time limit of the covariance a\ (t) is entirely 
determined by the asymptotic behavior of Ca(t). This 
function depends explicitely on the observable A(p). It 
is worthwhile to notice that ergodicity will therefore in 
general depend on the dynamical variable A{p) under 
consideration. In the following, we take A(p) = p n 
and evaluate the long-time behavior of Ca(t) follow- 
ing Ref. We find Ca(t) ~ with an expo- 
nent ll = (1 - (2n + l)D)/2D. As a result, the co- 
variance IjlUI) will converge to zero as t — > oo, only if 
D < D n = l/(2n + 1). We thus obtain the important 
result that there is a noise threshold D n above which er- 
godicity is broken. As already mentioned, this threshold 
depends explicitely on the parameter n, that is, on the 
quantity A{p): the smaller the power n, the larger the 
value of D n . On the other hand, the moments of the sta- 
tionary momentum distribution, (p m ) — J dpp m W s (p), 
are finite for D < D' m = l/(m+l). We can therefore con- 
clude that there exists a direct relationship between the 
loss of ergodicity in the system for the variable A(p) = p n 
and the divergence of the 2nth moment of the stationary 
momentum distribution. 

Let us look in more detail at the value n = 0, which 
corresponds to the largest noise threshold D n= o = 1. We 
first mention that for D > D n= o, the Fokker-Planck equa- 
tion doesn't have a normalizable stationary solution any- 
more and the system is obviously nonergodic. Further, 
for n — 0, the function Ca{t) can be rewritten in terms 
of the conditional probability density P 2 (p 2 , T~\pi, 0) as 



Ca{t) 



dpi W s {pi 



dp 2 (p 2 {p 2 ,r\p u Q)-W s {p 2 )) . (12) 



This result is of special interest. Equation (|12|) shows 
that, unlike for the Ornstcin Uhlcnbcck process, the 
steady state is here reached in a nonexponential way. 
Such an asymptotic power-law decay is usually associ- 
ated with non-Markovian processes and is rather sur- 
prising for a stationary Markov process as described by 
Eq. (JTJ. The origin of this algebraic behavior is of course 
rooted in the nonlinear drift coefficient @. 
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FIG. 1: Probability to be in the low-momentum region (con- 
tinuous line) and in the high-momentum region (dashed line) , 
in the long-time limit, as a function of the noise strength D. 

A closely related quantity is the first-passage time dis- 
tribution. The first-passage time is defined as the time 
at which the momentum of the system first exits a cer- 
tain momentum interval, given that it was originally in 
that interval. The first-passage time problem for the 
Fokker-Planck equation can again be treated using 
standard techniques of stochastic calculus (see for ex- 
ample Ref. |2lJ). In the high-momentum region, the 
Laplace transform g 2 (s) of the first-passage time distri- 
bution obeys the following backward equation, 



D 



d 2 g 2 1 dg 2 
dp 2 p dp 



- sg 2 = . 



(13) 



Solving Eq. i|13|l with the boundary conditions g(l) = 1 
and g(oo) = 0, we obtain 



92 (s) 



K v (py/I/D) 



P 



(14) 



where K u (p) is the modified Bessel function of the sec- 
ond kind of order v. It follows from Eq. (|14p. that 
to leading order, g 2 (s) ~ s v as s — > 0. The first- 
passage time distribution in region 2 is thus also a power- 
law tail distribution and it asymptotically behaves as 
g 2 (t) - i"T with an exponent 7 = (3D + 1)/{2D). The 
corresponding moments (t n ) = J dtt n g 2 (t) converge if 
D < D[[ = l/(2n — 1). The first-passage time distri- 
bution g\{t) in region 1 can be computed along similar 
lines and the associated moments can be shown to be all 
finite. Figure Q shows the probability to be in region 
1 and 2, in the limit of long times, as a function of the 
noise intensity D. We observe that for small D, atoms 
are mostly located in the low-momentum region, where 
they experience the linear part of the drift. On the other 
hand, for D close to D n= o, atoms get localized in the 
high-momentum region, where the drift asymptotically 
decays as —l/p. Remarquably, atoms in this system can 
thus be brought in a high-energy state through the sole 
action of the noise. We note that this problem exhibits 
an interesting analogy with subrecoil laser cooling, where 
atoms accumulate in a low-energy state (in a so-called 
dark state) , with infinite mean trapping time |22| . To our 
knowledge, loss of ergodicity in systems with divergent 
trapping times has been first discussed in the context of 
spin-glasses p3|. 
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We can now formulate the main result of the paper: 
Ensemble average and time average of the dynamical 
variable A(p) = p n stop being equal in the infinite-time 
limit — ergodicity of the system is accordingly broken — 
when the 2nth moment of the stationary momentum dis- 
tribution and the (n + l)th moment of the first passage 
time distribution in the high-momentum region become 
infinite. An unambiguous correspondence between the 
nonergodic properties of a system described power-law 
distributions and the divergence of their respective mo- 
ments is therefore demonstrated. This confirms and ex- 
tends the findings reported in Refs. H U |H |23|. We 
stress that the above result is not restricted to transport 
in an optical lattice and that our analysis applies to a 
whole class of systems described by an equation of the 
form Q with a drift coefficient decaying asymptotically 
as — 1/p. More generally, we conjecture that the result 
also holds true for other power-law systems, even if they 
are not described by a simple kinetic equation like the 
ordinary Fokker-Planck equation (JTJ) (llaEll^. 

Further insight can be gained by considering a dis- 
cretized form of the Fokker-Planck equation Q . In do- 
ing so, we shall obtain a generalization of the Ehrenfest 
urn model, which has played an important role in clar- 
ifying the foundations of statistical mechanics |24j. We 
write p = jAp and t = I At and find that the probability 
w(j, I) = W(jAp, I At) satisfies the difference equation, 

cj(j,l + l) = a(j-l)u(j-l,l) + b(j + l)u(j + l,l) . (15) 

In the limit (Ap, At)^0, Eq. 115JI reduces to the contin- 



uous equation <£Q) with K(p) — Ap/ At [a(p) — b(p)} and 
D= (Ap) 2 /(2At)[a(p) + b(p)}. The transition probabil- 
ities a(j) and b(j) in (|15|l are explicitely given by 



R(l+j 2 )-j 
2R(l + p) ■ 



b(j) 



2R(l + p) 



(16) 



We recall that the Ehrenfest model consists of two urns 
containing a total of 2R balls. At regular time intervals, 
At, a ball is randomly chosen and moved to the other 
urn; w(j, I) is then the probability of having R + j balls 
in the first urn at time I. In the standard version of the 
model, the j 2 terms in Eq. I|16l) are absent and Eq. 
reduces to the Ornstein-Uhlenbeck process with linear 
drift. We note that in the continuous description, the 
momentum p corresponds to the number of excess balls 
j in the first urn. So, for example, ergodicity breaking 
for n — 1 (D > 1/3) occurs when the fluctuations of the 
number of balls in each urn diverge, the average number 
of balls being still finite and equal. Moreover, when D 
approaches one, all the balls preferentially occupy the 
same urn, therefore acting as Maxwell's demon. 

In conclusion, we have investigated anomalous trans- 
port in an optical lattice from the point of view of sta- 
tistical mechanics and established an explicit correspon- 
dence between ergodicity breaking and the divergence of 
the moments of the power-law tail distributions describ- 
ing the behavior of the system, both in momentum space 
and in time. 
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